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We analyze observational constraints from TeV astrophysics on Lorentz violating nonlinear dis¬ 
persion for photons and electrons without assuming any a priori equality between the photon and 
electron parameters. The constraints arise from thresholds for vacuum Cerenkov radiation, photon 
decay and photo-production of electron-positron pairs. We show that the parameter plane for cubic 
momentum terms in the dispersion relations is constrained to an order unity region in Planck units. 
We find that the threshold configuration can occur with an asymmetric distribution of momentum 
for pair creation, and with a hard photon for vacuum Cerenkov radiation. 
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There are several reasons to suspect that Lorentz in¬ 
variance may be only a low energy symmetry. This possi¬ 
bility is suggested by the ultraviolet divergences of local 
quantum field theory, as well as by tentative results in 
various approaches to quantum gravity and string the¬ 
ory 0 11,11- Moreover, Lorentz symmetry can only 
ever be verified up to some finite observationally accessi¬ 
ble velocity, which leaves untested an infinite volume of 
the supposed symmetry group. 

The possibility of Lorentz violation can be studied, 
without a particular fundamental theory in hand, by con¬ 
sidering its manifestation in dispersion relations for par¬ 
ticles. If rotational invariance is preserved, it is natural 
to assume that deviations from the Lorentz invariant dis¬ 
persion relation (p) = -I- can be characterized 

at low energies by an expansion with integral powers of 
momentum, = m? -\- o-nP^, where the a„ 

are coefficients with mass dimension 2 — n which might 
be positive or negative. [Throughout this letter p de¬ 
notes the absolute value of the 3-momentum vector p, 
and we use units with the low energy speed of light in 
vacuum equal to unity.] Different approaches to quantum 
gravity suggest different leading order Lorentz violating 
terms. The terms with n < 4 have mostly been consid¬ 
ered so far. Observations limit the coefficients ai _2 to be 
extremely small (see e.g. §0,1 and references therein). 
In this letter we shall assume they are precisely zero. 

The cubic and higher order coefficients have negative 
mass dimension, so if the Lorentz violation descends from 
the Planck scale Mp = {hc^ ~ 1.22 • 10^® GeV, a„ 

would be expected to be of order Mp~^. In this case 
the Lorentz violation is naturally suppressed, and the 
lowest order term would dominate all the higher ones 
as long as the momentum is less than Mp, hence we 
shall restrict to a single Lorentz violating term. We thus 
consider the constraints that high energy observations 
impose on dispersion relations of the form 

El=pl + ml + riaPa/M^-‘^, ( 1 ) 

where a labels different fields and n > 3. We have intro¬ 


duced the energy scale M = 10^® GeV Mp so that the 
coefficients pa are dimensionless. If the Lorentz violation 
comes from quantum gravity effects, one would expect 
Pa to be within a few orders of magnitude of unity. In 
the absence of a fundamental theory one has no reason to 
expect any particular relation between the coefficients pa 
for different particles, except perhaps that they should all 
be of the same order of magnitude. Since the dispersion 
relation ( 0 ) is not Lorentz invariant, it can only hold in 
one reference frame. We assume along with all previous 
authors that this frame coincides with that of the cosmic 
microwave background. The velocity of the earth relative 
to this frame is negligible for our purposes. 

Observational consequences of the p term in (|^) may 
seem out of reach because of the Planck scale suppres¬ 
sion. However this is not so. Dispersion relations like 
(Q) produce kinematic relations from energy-momentum 
conservation that differ from the usual Lorentz invariant 
case. As a result reactions can take place that are nor¬ 
mally forbidden, and thresholds for reactions are modi¬ 
fied. One can expect deviations from standard threshold 
kinematics when the the last two terms of (0) are of com¬ 
parable magnitude. Assuming p is of order unity this 
yields the condition pdev ~ (m^M"-^)!/"^ which is ~ 
(rua/me)^/^ X 10 TeV for n = 3 and ~ (ma/we)^/^ x lO'^ 
TeV for n = 4, where rrie is the electron mass. Al¬ 
though these energies are currently not achievable in par¬ 
ticle accelerators (except in the case of massive neutri¬ 
nos which however are too weakly coupled to provide 
constraints) they are in the range of current astrophys- 
ical observations. In fact, it has been suggested by sev¬ 
eral authors ||, ^ |ll[ |^, |l^, |4[ |l^ (see also Q 
and references therein) that we may already be observ¬ 
ing deviations from Lorentz invariance via the possibly 
missing Greisen, Zatsepin and Kuzmin (GZK) cut-off 
on cosmic ray protons with ultra high energy greater 
than 7 x 10^® eV [|^, and the possible overabundance 
of gamma rays above 10 TeV from the blazar system 
Markarian 501 Here we shall mostly not con¬ 

sider the constraints imposed by asking Lorentz violation 
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to explain these puzzles. Instead we restrict our atten¬ 
tion to constraints imposed by consistency with known 
phenomena (or lack thereof). 

Observational constraints'. Several studies of obser¬ 
vational limits on Lorentz violating dispersion relations 
have already been carried out | 0 , ||, ||, |^, |^, |^, |^, [l^ , 
with various different assumptions about the 
coefficients. Our study focuses on purely QED interac¬ 
tions involving just photons and electrons. We assume 
n = 3, since the n = 4 terms are suppressed by another 
inverse power of M. Unlike other studies, no a priori 
relation between the coefficients r].y and rje is assumed, 
and we combine all the different constraints in order to 
determine the allowed region in the parameter plane. To 
eliminate the subscript a we introduce ^ := 77 -,,, rj := rje, 
and m := me- 


The modified dispersion relations for photons and elec¬ 
trons in general allow two processes that are normally 
kinematically forbidden: vacuum Cerenkov radiation, 
e~ —> e“ 7 , and photon decay, 7 —> e+e”. In addition 
the threshold for photon annihilation, 77 —> e“'"e“, is 
shifted. The vacuum Cerenkov process is extremely ef¬ 
ficient, leading to an energy loss rate that goes like 
well above threshold. Similarly the photon decay rate 
goes like E. Thus any electron or photon known to prop¬ 
agate must lie below the corresponding threshold. 


We consider constraints that follow from three consid¬ 
erations: (i) Electrons of energy ~ 100 TeV are belived 
to produce observed X-ray synchrotron radiation com¬ 
ing from supernova remnants [^, and to also produce 
multi-TeV photons by inverse-Compton scattering with 
these X-rays |^. Assuming these electrons are ac¬ 


tually present, vacuum Cerenkov radiation must not oc¬ 
cur up to that energy (ii) Gamma rays up to 50 
TeV arrive on earth from the Crab nebula |^, so pho¬ 
ton decay does not occur up to this energy, (iii) Cosmic 
gamma rays are believed to be absorbed in a manner 
consistent with photon annihilation off the infrared (IR) 
background with the standard threshold |Q. Observa¬ 
tion (iii) is not model independent, so the corresponding 
constraint is tentative and subject to future verification. 


Modified kinematics'. The processes e~ —> e “7 and 
7 —> e+e” correspond to the basic QED vertex, but are 
normally forbidden by energy-momentum conservation 
together with the standard dispersion relations. When 
the latter are modified, these processes can be allowed. 


^ The competing energy loss by synchrotron radiation is irrelevant 
for this constraint. The rate of energy loss from a particle of 
energy E due to the vacuum Cerenkov effect goes like —e^E^, 
while that from synchrotron emission goes like —e'^B^E^/m'^ 
(using units where c = 7i = 1). For a magnetic field of about 
one micro Gauss (as those involved in supernova remnants) the 
synchrotron emission rate is 40 orders of magnitude smaller than 
the vacuum Cerenkov rate. 


To see this, let us denote the photon 4-momentum by 
ki = (wfcjk), and the electron and positron 4-momenta 
by Pa = {Ep,p) and = {Eq,q). Eor the two reactions 
energy-momentum conservation then implies P 4 = k 4 + q 4 
and k 4 — P 4 + 54 respectively. In both cases, we have 
{p 4 — k 4 )‘^ = ql, where the superscript “ 2 ” indicates the 
Minkowski squared norm. Using the Lorentz breaking 
dispersion relation Eq. (Q) this becomes 

fk^ + pp^ — pq^ = 2M {EpUik — pk cos 0 ), (2) 


where 0 is the angle between p and k. In the standard 
case the coefficients ^ and p are zero and the r.h.s. of 
Eq. (^) is always positive, hence there is no solution. It 
is clear that non-zero ^ and p can change this conclusion 
and allow these processes. 

To derive the observational constraints one needs to 
determine the threshold for each process, i.e. the lowest 
energy for which the process occurs. Assuming mono¬ 
tonicity of all the dispersion relations (for the relevant 
momenta ^ M) one can show that all thresholds 
for processes with two particle final states occur when 
the final momenta are parallel. Moreover for two parti¬ 
cle initial states the incoming momenta are antiparallel. 
This implies that at a threshold 0 = 0 in Eq. (||) and that 
in the corresponding formula for the photon annihilation 
we shall consider antiparallel photons and parallel lep¬ 
tons. These geometries have been assumed in previous 
works but to our knowledge they were not shown to be 
necessary. In fact they are not necessary if the disper¬ 
sion relations are not monotonic. Details concerning the 
determination of the thresholds are reported in . 

Vacuum Cerenkov radiation: We find that an electron 
can emit Cerenkov radiation in the vacuum if 77 > 0 or 
if ?7 < 0 and ^ < p. Depending on the values of the 
parameters, the threshold configuration can occur with a 
zero-energy photon or with a finite energy photon. These 
two cases correspond to the following two threshold rela¬ 
tions, respectively: 


Pth 


Pth 



for 77 > 0 and ^ > —3p, (3) 


AmfM (^ -I- 77 ) 

(C - vf 


for ^ < -377 < 0, 

(4) 

or f < p <0. 


The reaction is not allowed in the region where f > p 
and 77 < 0. Note that if ^ = 77 only the solution (^ yields 
a finite threshold. 

Electrons of energy ~ 100 TeV are indirectly observed 
via X-ray synchrotron radiation coming from supernova 
remnants P2[ . Thus for example in the region of the 
parameter plane where (H) holds we obtain the constraint 
77 < m?M/2p\^ ^ 10 “^. 

Photon decay: A photon can spontaneously decay into 
an electron-positron pair provided ^ is sufficiently great 








3 


for any given rj. Contrary to Lorentz-invariant kinematics 
of pair creation thresholds, we find that the two particles 
of the pair do not always have equal momenta. Photon 
decay is allowed above a broken line in the rj-^ plane 
given by ^ = 77/2 in the quadrant 77 > 0 and by ^ = 77 
in the quadrant 77 < 0. Above this line, the threshold 


is given by 







for C > 0, 

(5) 

"-th 



^th — 

—Sm^Mr] 

1/3 

for 77 < ^ < 0. 

(6) 

The first relation (||) arises 

when the electron 

and 


positron momenta are equal at threshold. The second 
relation (|^) applies in the case of asymmetric distribu¬ 
tion of momenta. Note that if ^ = 77 , the asymmetric 
threshold disappears, leaving just the symmetric one. 

The constraint we impose is that the threshold is above 
50 TeV, the highest energy of observed gamma rays from 
the Crab nebula |^. The strength of the constraint is 
determined by the smallness of the quantity 
For fcniax = 50 TeV one gets « 0.02. 

Photon annihilation: The standard threshold for a 
gamma ray to annihilate with an IR background photon 
of energy e is = Tn^/e. In the presence of dispersion the 
threshold relations take approximately the same form as 
for photon decay, equations (H,||), with the replacement 
^ ^ , where ^ (Here we have used the 

fact that e is much smaller than any other scale in the 
problem.) However, now these relations correspond re¬ 
spectively to cubic and quartic polynomial equations for 
fcth (since is itself a function of and the condition 
that determines whether the threshold is given by the 
symmetric (^) or asymmetric (^) relation is more com¬ 
plicated. The detailed analysis can be found in Here 
we merely state the result. Rather than fixing 77 , ^ and 
e and solving the relations for fcth, we fix e and fcth and 
use the threshold relations to solve for ^ as a function of 
77 . When fcth < 1.5/cs the symmetric threshold applies for 

> 0 and the asymmetric one applies for 77 < ^' < 0 . 
When fcth > 1.5fcs there is no symmetric threshold, and 
the asymmetric one applies only for ^ below the symmet¬ 
ric fcth = 1.5fcs line, ^ = 77/2 — {16/27){e^M/m'^). In the 
case ^ = rj the threshold configuration is never asymmet¬ 
ric 

For the observational consequences it is important to 
recognize that the threshold shifts are much more signifi¬ 
cant at higher energies than at lower energies. To exhibit 
this dependence, it is simplest to fix a gamma ray energy 
k and to solve for the corresponding soft photon thresh¬ 
old energy eth ■ Taking the ratio with the usual threshold 
eth, 0 i we find a dependence on k at least as strong as 
Introducing fcio := fc/(10TeV), we have 


eth 

eth,o 


(77^^ 3 

^ + on ''lO 


for e' > 0, (7) 


“ = ^10 + \/-T^w for <0- (8) 

eth,o V u 

High energy TeV gamma rays from the blazars Markar- 
ian 421 and Markarian 501 have been detected out to 
17 TeV and 24 TeV respectively |^, Although the 
sources are not well understood, and the intergalactic IR 
background is also not fully known, detailed modeling 
shows that the data are consistent with some absorp¬ 
tion by photon annihilation off the IR background (see 
e.g. H, H and references therein). However, while 
the inferred source spectrum for Markarian 501 is consis¬ 
tent with expectations for energies less than around 10 
TeV, above this energy there have been claims ||l^, |l^ 
that far more photons than expected are detected. Nev¬ 
ertheless, recent analysis based on a more detailed recon¬ 
struction of the IR background do not seem to corrobo¬ 
rate this point of view ||26| . 

Due to these uncertainties sharp constraints from pho¬ 
ton annihilation are currently precluded. Instead, we 
just determine the range of parameters 77 for which the 
threshold fcth lies between 10 TeV and 20 TeV for an IR 
photon of energy 0.025 eV with which a 10 TeV pho¬ 
ton would normally be at threshold. Based on current 
observations it seems unlikely that the threshold could 
lie far outside this range. (It has previously being pro¬ 
posed that raising this threshold by a factor of two 
could explain the potential overabundance of photons 
over 10 TeV.) Given the strong energy dependence of 
the threshold shift in equations (^ and (^) this thresh¬ 
old raising would not be obviously in disagreement with 
current observations below 10 TeV. 

Combined constraints: Putting together all the con¬ 
straints and potential constraints we obtain the allowed 
region in the 77-^ plane (see Figure The photon decay 
and Cerenkov constraints exclude the horizontally and 
vertically shaded regions, respectively. The allowed re¬ 
gion lies in the lower left quadrant, except for an exceed¬ 
ingly small sliver near the origin with 0 < 77 < 10“^ and 
a small triangular region (—0.16 ^77<0, 0 <^;$ 0.08) 
in the upper left quadrant. The range of the photon 
annihilation threshold previously discussed falls between 
the two roughly parallel diagonal lines. The upper diag¬ 
onal line corresponds to the standard threshold kg = 10 
TeV and the lower diagonal line to not more than twice 
that threshold. If future observations of the blazar fluxes 
and the IR background confirm agreement with stan¬ 
dard Lorentz invariant kinematics, the region allowed by 
the photon annihilation constraint will be squeezed to¬ 
ward the upper line (fcth = kg). This would close off all 
the available parameter space except for a region much 
smaller than unity around the Lorentz-invariant values 
$ = 0 = 77. 

Conclusions: We have shown that astrophysical ob¬ 
servations put strong constraints on the possibility of 
Lorentz-violating Planck scale cubic modifications to 
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the electron and photon dispersion relations. The con¬ 
straints arise due to the effect these modifications have on 
thresholds for various reactions. We have also seen that 
the threshold configurations with a final state electron- 
positron pair sometimes involve unequal momenta for the 
pair, unlike what occurs for all Lorentz-invariant decays. 
This can happen if ^ ^ 77 and ^,77 < 0. 

The allowed region in the 77 — ^ plane includes ^ = rj = 
— 1 , which has been a focus of previous work 
The negative quadrant has most of the allowed parameter 
range. Note that in this quadrant all group velocities are 
less than the low energy speed of light. 

To further constrain the cubic case will require new 
observations. Finding higher energy electrons would not 
help much, while finding higher energy undecayed pho¬ 
tons would squeeze the allowed region onto the line ^ = 77 . 
To shrink the allowed segment of this line using the re¬ 
actions we have considered would require observations 
confirming the usual threshold for photon annihilation 
to higher precision. 

Perhaps other processes could be used as well. One 
might have hoped that observations comparing the time 
of flight of photons of different frequencies from distant 
sources such as gamma ray bursts and active galactic 
nuclei would help constrain the absolute value of ^ (see 
e.g. I, il H). Unfortunately current observations 
just yield |^| < 122 for 71 = 3. This is an interest¬ 
ing constraint but it is not competitive with the other 
ones already considered here. (However the forthcom¬ 
ing Gamma Ray Large Area Space Telescope (GLAST) 
mission may provide more stringent constraints of this 
type [Q.) Another idea is to exploit the fact that the re¬ 
action 7 —> 87 is kinematically allowed with finite phase 
space and nonzero amplitude in the presence of modified 
dispersion, unlike in the standard case. This photon de¬ 
cay channel occurs at all energies if ^ > 0 , i.e. it has no 
threshold, so it might be thought to provide a very pow¬ 
erful constraint on positive Unfortunately, however, 
the amplitude for this reaction is far too small to provide 
any useful constraint [ p^ . 

It is interesting to consider the case of the possibly 
missing GZK cutoff . If the cutoff is really missing, it 
has been proposed to explain this using Lorentz violating 
dispersion &§■ The relevant protons are at such a high 
energy — over 10^® eV — that it takes only tiny Lorentz 
violating parameters rja in (|^) to increase the threshold 
by an amount of order unity or more. In particular, if one 
assumes all coefficients rja are equal, this only requires 77 
negative and I77I > /p'^ ~ ^Q-38+9n^ Ppj. ^ = 3 

this is 10“^^, and for 77 = 4 it is still only 10“^. Thus for 
both the 77 = 3 and 77 = 4 cases only very small values of 
77 are needed to dramatically modify the GZK cutoff, so 
a shifted cutoff could be explained by Lorentz-violating 
constants with our constraints. However recent data |3^, 
^ strongly support the existence of the GZK cutoff at its 
expected (Lorentz invariant) value. If this is confirmed. 


the above analysis shows that the GZK reaction provides 
very good constraints for modifications up to 77 = 4 p8|. 
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FIG. 1: Combined constraints on the dimensionless photon 
and electron parameters for the case 77 = 3 ( 77 , ^ = 1 corre¬ 
sponds to the dimensionful coefficient M~^ = 10“^® GeV“^ 
of the 77 = 3 term in Eq. ©)• The regions excluded by the 
photon decay and Cerenkov constraints are lined horizontally 
in blue and vertically in red respectively. The region between 
the two diagonal green lines corresponds to a threshold be¬ 
tween one and two times the standard threshold (which is 
10 TeV for photon annihilation with an IR photon of energy 
0.025 eV). The upper green line corresponds to the unmodified 
threshold. The shaded patch is the part of the allowed region 
that falls between these photon annihilation thresholds. The 
dashed line is ^ = 77 . 
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